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' Abstract. 

i Hydrodynamic description requires a local thermodynamic equilibrium of the system under 

study but an approximate hydrodynamic behaviour is already manifested when a momentum 
CN| ' distribution of liquid components is not of equilibrium form but merely isotropic. While the 

process of equilibration is relatively slow, the parton system becomes isotropic rather fast due to 
the plasma instabilities. Azimuthal fluctuations observed in relativistic heavy-ion collisions are 
argued to distinguish between a fully equilibrated and only isotropic parton system produced 
in the collision early stage. 



1. Introduction 

A matter created in relativistic heavy-ion collisions manifests a strongly collective hydrodynamic 
behaviour [J, particularly evident in studies of the so-called elliptic flow [2J |3J HJ 
Hydrodynamic description requires, strictly speaking, a local thermal equilibrium and 
experimental data on the particle spectra and the elliptic flow suggest, when analysed within 
the hydrodynamic model, that an equilibration time of the parton 2 system produced at the 
collision early stage is as short as 0.6 fm/c [5]. Such a fast equilibration can be explained 
assuming that the quark-gluon plasma is strongly coupled 7 . However, high-energy density 
in the collision early stage, when the elliptic flow is generated jHj, allows one to believe that 
the plasma is then weakly coupled due to the asymptotic freedom. Calculations, which assume 
that the parton-parton collisions are responsible for the equilibration of the weakly interacting 
plasma, provide an equilibration time of at least 2.6 fm/c (Sj. To thermalize the system one 
needs either a few hard collisions of the momentum transfer of order of the characteristic parton 
momentum 3 , which we denote here as T (as the temperature of equilibrium system), or many 
collisions of smaller transfer. As discussed in e.g. ^Oj; the inverse time scale of the collisional 
equilibration is of order g i \n(\/g)T where g is the QCD coupling constant. However, it has 
been argued that the equilibration is speeded up by instabilities generated in an anisotropic 
quark-gluon plasma as growth of the unstable modes is associated with the system's 
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2 The term 'parton' is used to denote a fermionic (quark) or bosonic (gluon) excitation of the quark-gluon plasma. 

3 Although I consider anisotropic systems, the characteristic momentum in all directions is assumed to be of the 
same order. 



isotropization. The characteristic inverse time of instability development is roughly of order 
gT for a sufficiently anisotropic momentum distribution I14 ( IT5 | I16j. Thus, the instabilities 
are much 'faster' than the hard collisions in the weak coupling regime. Very recent classical 
simulation indeed shows effectiveness of the instabilities driven isotropization. 

The isotropization should be clearly distinguished from the equilibration process 
The instabilities driven isotropization is a mean-field reversible phenomenon which is not 
accompanied with entropy production. Therefore, the collisions, which are responsible for the 
dissipation, are needed to reach the equilibrium state of maximal entropy. The instabilities 
contribute to the equilibration indirectly, reducing relative parton momenta and increasing the 
collision rate. 

It has been recently observed that the hydrodynamic collective behaviour does not actually 
require local thermodynamic equilibrium but a merely isotropic momentum distribution of liquid 
components |13j . Thus, there is a question whether a quark-gluon plasma, which is equilibrated 
nearly immediately after its production as advocated in [7j, can be distinguished from the 
parton system which slowly evolves towards equilibrium being isotropized fast. I argue here 
that measurements of azimuthal fluctuations, which are generated at the early stage of heavy- 
ion collisions, can help to distinguish the two scenarios. 

In the first part of my talk I review the instabilities driven isotropization. I discuss how the 
unstable modes are initiated and what is a mechanism responsible for their growth. Dispersion 
relations of the unstable modes are considered, and it is explained why the development of 
instabilities is associated with the system's isotropization. In the second part of my talk I 
discuss the azimuthal fluctuations, arguing that the fluctuations generated in the non-equilibrium 
isotropic system are much larger than those in the fully equilibrated plasma. Two possible 
measurements are proposed. 

2. Instabilities driven isotropization 

Temporal evolution of the electron-ion plasma is plagued by a large variety of instabilities. Those 
caused by coordinate space inhomogeneities, in particular by the system's boundaries, are usually 
called the hydrodynamic instabilities while those due to non-equilibrium momentum distribution 
of plasma particles the kinetic instabilities. Hardly anything is known about hydrodynamic 
instabilities of the quark-gluon plasma, and I will not speculate about their possible role in 
the system's dynamics. The kinetic instabilities are initiated either by the charge or current 
fluctuations. In the first case, the electric field (E) is longitudinal (E || k, where k is the wave 
vector), while in the second case the field is transverse (E _L k). For this reason, the kinetic 
instabilities caused by the charge fluctuations are usually called longitudinal while those caused 
the current fluctuations transverse. Since the electric field plays a crucial role in the longitudinal 
mode generation, the longitudinal instabilities are also called electric while the transverse ones 
magnetic. 

In the non-relativistic plasma the electric instabilities are usually much more important 
than the magnetic ones as the magnetic effects are suppressed by the factor v 2 /c 2 where 
v is the particle's velocity. In the relativistic plasma both types of similar strength. The 
electric instabilities occur when the momentum distribution has more than one maximum while 
a sufficient condition for the magnetic instabilities is, as discussed below, anisotropy of the 
momentum distribution. For this reason, the magnetic unstable mode, which is also called 
Weibel or filamentation instability JB]> was argued long ago to be relevant for equilibration of 
the quark-gluon plasma produced in relativistic heavy-ion collisions In the remaing part of 
the section, I am going to explain in detail why the filamentation is relevant and how it speeds 
up the process of plasma thermalization. 



2. 1 . Seeds of the filamentation 

Let me first discuss how the unstable transverse modes are initiated. For this purpose I consider 
a parton system which is homogeneous but the parton momentum distribution is, in general, 
not of the equilibrium form, it is not isotropic. The system is on average locally colourless but 
colour fluctuations are possible. Therefore, (j£(x)) = where j£(x) is a local colour four-current 
in the adjoint representation of SU(3) gauge group with \x = 0, 1, 2, 3 and a = 1, 2, 3, . . . , 8 being 
the Lorentz and colour index, respectively; x = (t, x) denotes a four-position in the coordinate 
space. 

As discussed in detail in |19j . the current correlator for a classical system of non- interacting 
quarks and gluons is 

M%fax) d ^ f ^(t 1> x 1 )j?(< 4> x a )> = i^^y"^^/(p)^)(x-vt), (1) 

where (t,x) = fa — t±,X2 — xi) and the effective parton distribution function /(p) equals 
n(p)+n(p)+6n 5 (p) with n(p), n(p) and n g (p) giving the average colourless distribution function 
of quarks Q^(x, p) = 5*%(p), antiquarks Q^(x,p) = 5 t: >n(p), and gluons G ab (x,p) = 5 ab n g (p). 
We note that the distribution function of (anti-)quarks belongs to the fundamental representation 
of the SU(3) gauge while that of gluons to the adjoint representation. Therefore, i, j = 1,2,3 
and a, b = 1, 2, 8. 

Due to the average space-time homogeneity, the correlation tensor Q depends only on the 
difference fa — ti,X2 — xi). The space-time points fa,Xi) and fa,X2) are correlated in the 
system of non-interacting particles if a particle travels from fa,xi) to (t2) x 2)- For this reason 
the delta <5( 3 )(x — vt) is present in the formula (jTJ). The momentum integral of the distribution 
function simply represents the summation over particles. The fluctuation spectrum is found as 
a Fourier transform of the tensor Q i- e. 

M£>, k) = l<? <?*J ^ ^ /(p) 2n5(u - kv) . (2) 

To compute the fluctuation spectrum, the parton momentum distribution has to be specified. 
Such calculations with two forms of the anisotropic momentum distribution are presented in 
|19j . Here I only qualitatively discuss Eqs. (1ll2f) . I assume that the momentum distribution is 
elongated in, say, the z direction. Then, Eqs. (1 1 121) clearly show that the corelator M zz is larger 
than M xx or M yy . It also clear that M zz is the largest when the wave vector k is along the 
direction of the momentum deficit. Then, the delta function 5{uo — kv) does not much constrain 
the integral in Eq. (J2J). Since the momentum distribution is elongated in the z direction, the 
current fluctuations are the largest when the wave vector k is the x—y plane. Thus, I conclude 
that some fluctuations in the anisotropic system are large, much larger than in the isotropic one 
and that anisotropic system has a natural tendency to split into the current filaments parallel to 
the direction of the momentum surplus. These currents are seeds of the filamentation instability. 

2.2. Mechanism of filamentation 

Let me now explain in terms of elementary physics why the fluctuating currents, which flow in the 
direction of the momentum surplus, can grow in time. To simplify the discussion, which follows 
|19j . I consider an electromagnetic anisotropic system. The form of the fluctuating current is 
chosen to be 

j(x) = j e z cos(k x x) , (3) 

where e z is the unit vector in the z direction. As seen in Eq. (j3J), there are current filaments 
of the thickness ir/\k x \ with the current flowing in the opposite directions in the neighbouring 
filaments. 




The magnetic field generated by the current Q is given as 

B(x) = 77- e y sm(k x x) , 

&x 

and the Lorentz force acting on the partons, which fly along the z direction, equals 

j 

F(x) = ?vx B(x) = -q v z — e x sin(k x x) , 

&x 

where q is the electric charge. One observes, see Fig. 1, that the force distributes the partons in 
such a way that those, which positively contribute to the current in a given filament, are focused 
in the filament centre while those, which negatively contribute, are moved to the neighbouring 
one. Thus, the initial current is growing. 



2.3. Dispersion equation 

The Fourier transformed chromodynamic field A^{k) satisfies the equation of motion as 

fcV - W - U^(k)] A v {k) = , (4) 

where IP" (A;) is the polarization tensor or gluon self-energy which is discussed later on. A 
general plasmon dispersion equation is of the form 



det 







(5) 



Equivalently, the dispersion relations are given by the positions of poles of the effective gluon 
propagator. Due to the transversality of 11^ (k) (k fJi H^ L ' (k) = k u Il^ u (k) = 0) not all components 
ofU^ u (k) are independent from each other, and consequently the dispersion equation which 
involves a determinant of 4 x 4 matrix, can be simplified to the determinant of 3 x 3 matrix. For 
this purpose I introduce the colour permittivity tensor e lm (k) where the indices l,m,n = 1,2,3 
label three-vector and tensor components. Because of the relation 



e lm (k)E l {k)E m (k) = W(k)A^k)A v (k) , 



where E is the chromoelectric vector, the permittivity can be expressed through the polarization 
tensor as 

JLm(i„\ j^lm i 1 Trim/ 



UJ 1 



Then, the dispersion equation gets the form 

det 



k 2 6 lm 



k l k m -uj 2 e lm (k) 



0. 



(6) 



The relationship between Eq. © and Eq. © is most easily seen in the Coulomb gauge when 
A = and k • A(/c) = 0. Then, E = iujA. and Eq. (JIJ) is immediately transformed into an 
equation of motion of E(/c) which further provides the dispersion equation ©. 

The dynamical information is contained in the polarization tensor Ii^ v {k) or, equivalently, 
in the permittivity tensor e lm (k) which can be derived either within the transport theory or 
diagrammatically j^Oj. The result is 



V,k) 



g 2 
+ — 

2uj 
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df( P ) 



(2vr) 3 uj - kv + i0+ dp 1 



l„,m 



k V 



UJ ' UJ 



(7) 



Since H fiu (k) and e lm (k) are unit matrices in the colour space, the colour indices are suppressed 
here. 

Substituting the permittivity ((7| into Eq. ©, one fully specifies the dispersion equation © 
which provides a spectrum of quasi-particle bosonic excitations. A solution w(k) of Eq. (jEJ) is 
called stable when Imu; < and unstable when Imu; > 0. In the first case the amplitude is 
constant or it exponentially decreases in time while in the second one there is an exponential 
growth of the amplitude. In practice it appears difficult to find solutions of Eq. © because 
of rather complicated structure of the tensor 0. However, the problem simplifies as we are 
interested in specific modes which are expected to be unstable. Namely, we look for solutions 
corresponding to the fluctuating current in the direction of the momentum surplus and the wave 
vector perpendicular to it. 

As previously, the momentum distribution is assumed to be elongated in the z direction, 
and consequently the fluctuating current also flows in this direction. The magnetic field has a 
non-vanishing component along the y direction and the electric filed in the z direction. Finally, 
the wave vector is parallel to the axis x, see Fig. 1. We also assume that the momentum 
distribution obeys the mirror symmetry /(— p) = f(p), and then the permittivity tensor has 
only non-vanishing diagonal components. Taking into account all these conditions, one simplifies 
the dispersion equation |HJ) to the form 



H(oj) = k A x - uj 2 e zz (oj, k x ) = 



(8) 



where only one diagonal component of the dielectric tensor enters. 

It appears that an existence of unstable solutions of Eq. © can be proved without solving it. 
The so-called Penrose criterion |22j . which follows from analytic properties of the permittivity 
as a function of uj, states that the dispersion equation H(uj) = has unstable solutions if 
H(uj = 0) < 0. The Penrose criterion was applied to the equation (jHJ) in but a much more 
general discussion of the instability condition is presented in |16j . Not entering into details, 
there exist unstable modes if the momentum distribution averaged (with a proper weight) over 
momentum length is anisotropic. 

To solve the dispersion equation ©, the parton momentum distribution has to be 
specified. Several analytic (usually approximate) solution of the dispersion equation with various 
momentum distributions can be found in [111 1151 116j . An example of the numerical solution, 



which gives the unstable mode frequency in the full range of wave vectors is shown Fig. 2 taken 
from Jl]. The momentum distribution is of the form 



where p± = 
coupling is a, 



pI+pI- 



f(p) 



(Pt + a 



2 \3 



The mode is pure imaginary and jk = lmuj(k±^). The value of the 
= 0.3, aj_ = 0.3 GeV and the effective parton density is chosen to be 



6 fm .As seen, there is a finite interval of wave vectors for which the unstable modes exist. 
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Figure 2. The growth rate of the unstable 
mode as a function of the wave vector k = 
(Xl,0,0) for a ± = 0.3 GeV and 4 values 
of the parameter cry which controls system's 
anisotropy. The figure is taken from |14j . 
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Figure 3. Time evolution of the kinetic 
energy of particles (upper panel) and the 
energy stored in electric and magnetic 
fields (lower panel) in GeV/fm 3 for the 
U(l) and SU(2) gauge groups. The 
figure is taken from |17j . 



2.4- Growth of instabilities and abelianization 

A time evolution of a classical many-parton system interacting via classical chromodynamic 
field has been studied in ^7j- Numerical simulations have been performed effectively in 1 + 1 
dimensions as the chromodynamic potentials depend on t and x. The initial field amplitudes 
are assumed to obey Gaussian white noise and the initial parton momentum distribution is 

\/pl+Py 

/(P) ~ HPz) e p ha rd , (9) 

With Phard = 10 GeV. 

Fig. 3, which is taken from |17j . shows results of the simulation corresponding to a lattice 
of physical size L = 40 fm. As seen, the amount of energy of the fields, which is initially 
much smaller than the kinetic energy of all particles, grows exponentially and the magnetic 
contribution dominates. The simulation J2j indirectly confirms existence of the unstable 
magnetic modes in the system. 



Figure 4. The effective potential of the unstable magnetic mode as a function of two colour 
components of A a belonging to the SU(2) gauge group. The figure is taken from [2*T] . 

Unstable modes cannot grow to infinity and the question arises what is a mechanism 
responsible for stopping the instability growth. One suspects that non-Abelian non-linearities 
can play an important role here. An elegant qualitative argument [2^ suggests that the non- 
linearities do not stabilize the unstable modes because the system spontaneously chooses an 
Abelian configuration in the course of the instability development. Let me explain the idea. 

In the Coulomb gauge the effective potential of the unstable configuration has the form 



which is shown in Fig. 4 taken from The first term (with [i 2 > 0) is responsible for a very 
existence of the instability. The second term, which comes from the Yang-Mills lagrangian, is of 
pure non-Abelian nature. The term appears to be positive and thus it counteracts the instability 
growth. However, the non-Abelian term vanishes when the potential A a is effectively Abelian, 
and consequently, such a configuration corresponds to the steepest decrease of the effective 
potential. Thus, the system spontaneously abelianizes in the course of instability growth. 

The effect of abelianization has been indeed found in numerical simulations performed in 
the 1 + 1 dimensions |12 [ 1171 l2*T] . As an example, I show in Fig. 5 the result of fully classical 
simulation ^Jj. One observes in Fig. 5 taken from where 



that the field commutator measured by C decreases in time, in spite of the field growth quantified 

by (Arms- 

Very recent simulations performed in the 1 + 3 dimensions |231 I24| . which utilise a complete 
Hard Loop action for anisotropic systems |25| . show that the growth of unstable modes, which is 
initially exponential, becomes only linear at the later times. And the abelianization works 
only in the exponential period of instability development. However, it might well be that 
the abelianization becomes more efficient when the dynamical effects beyond the Hard Loop 
approximation are taken into account |26| . 

2.5. Isotropization 

When instabilites grow the systems becomes more isotropic because the Lorentz force acts on 
particle's momenta and the growing fields generate an extra momentum. 
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Figure 5. Temporal evolution of C and Figure 6. Temporal evolution of T xx and 
f/> rms measured in GeV. The figure is taken (T yy + T zz )/2. The figure is taken from [T7] . 
from (17j . 

To explain the mechanism I assume, as previously, that initially there is a momentum surplus 
in the z direction. The fluctuating current tends to flow in the z direction with the wave vector 
pointing in the x direction. Since the magnetic field has a y component, the Lorentz force, which 
acts on partons flying along the z axis, pushes the partons in the x direction where there is a 
momentum deficit. 

The effect of isotropization due to the action of the Lorentz force is nicely seen in the classical 
simulation ^7J. In Fig. 6, which is taken from ^7], there are shown diagonal components of the 
energy-momentum tensor 



The initial momentum distribution is given by Eq. Q, and consequently T xx = at t = 0. As 
seen in Fig. 6, T xx exponentially grows. 

The system isotropizes not only due to the effect of the Lorentz force but also due to the 
momentum carried by the growing field. As explained in detail in the momentum of the 
field is oriented along the wave vector which points in the direction of the momentum deficit. 
This effect has not been numerically studied yet but it is clear that the effect is comparable to 
that of Lorentz force only for suffuciently large field amplitudes. 

3. Azimuthal fluctuations 

In the first part of my talk I have argued that the quark-gluon plasma becomes isotropic fast 
due to the magnetic instabilities. And it has been recently observed ^3] that the system 
with isotropic momentum distribution manifests a hydrodynamic collective behaviour. The 
question arises whether such an approximate hydrodynamics can be distinguished from the real 
hydrodynamics describing a system which is in a local thermodynamic equilibrium. In the 
second part of my talk I propose to address the question by studying the azimuthal fluctuations. 

In relativistic heavy- ion collisions both at CERN SPS and BNL RHIC, one observes a sizable 
elliptic flow which is quantified by the second angular harmonics v<i of the azimuthal distribution 
of final state hadrons 12 EJ O2] • The phenomenon, which is sensitive to the collision early 
stage JHJ when the interaction zone is of the almond shape, is naturally explained within a 
hydrodynamics as a result of large density gradients [27J |2Hl E2 IHOj • Hydrodynamic description 
requires that the system under study is in a local thermodynamical equilibrium. However, 
an approximate hydrodynamic behaviour occurs, as argued in ^2], when the momentum 
distribution of liquid components is merely isotropic in the local rest frame. The point is 




that the structure of the ideal fluid energy- momentum tensor i.e. T^ v = (e + p)u^u v — pg^ u , 
where e, p and u M is the energy density, pressure and hydrodynamic velocity, respectively, 
holds for an arbitrary, though isotropic momentum distribution, e and p are then not the 
energy density and pressure but the moments of the distribution function which are equal the 
energy density and pressure in the equilibrium limit. Since the tensor T^ v obeys the continuity 
equation d^T^ v = 0, one gets an analogue of the Euler equation. However, due to the lack of 
thermodynamic equilibrium there is no entropy conservation and the equation of state is missing. 

Usually, non-equilibrium fluctuations are significantly smaller than the equilibrium 
fluctuations of the same quantity. A specific example of such an situation has been discussed 
in Sec. 2.1. Therefore, I expect that the fluctuations of V2 produced in the course of real 
hydrodynamic evolution are significantly smaller than those generated in the non-equilibrium 
quark-gluon plasma which is merely isotropic. It should be stressed here that the elliptic flow is 
generated in the collision early stage. Thus, I propose to carefully measure the fluctuations of 
t>2 as discussed in |31| . Since such a measurement is rather difficult, I also consider an integral 
measurement of azimuthal fluctuations proposed in [32] which can also help to distinguish the 
equilibrium from non-equilibrium fluctuations. 

3.1. Elliptic flow fluctuations 

In my discussion of V2 fluctuations I follow [JJ where the standard method |33l I34j to measure 
the elliptic flow was used. The method focuses on the angular distributions relative to direction 
of the impact parameter. The experimental procedure splits in two steps which should be as 
independent as possible. In the first step, one determines the impact parameter direction ipR, 
while in the second step, one constructs the distribution of azimuthal angle relative to ipn and 
one computes the Fourier coefficients. 

The one-particle distribution in a single event can be written as 



n=l 

Since the reaction plane is never reconstructed precisely and the real reaction plane angle ipR 
deviates from the estimated angle ipE, the n— th Fourier amplitude v n is determined as 



where R n = cos(n(tpR — ip E )) is the reaction plane resolution factor and 7TT denotes averaging 
over particles from a single event. 

Let me now consider an ensemble of events with every event representing a single nucleus- 
nucleus collision at fixed value (not direction) of the impact parameter. The angle ipR (and 
iPe) obviously varies form event to event. The question is how to detect the event-by-event 
fluctuations of the Fourier amplitudes v n . According to [31], the amplitude averaged over events 
is defined as 




(10) 



— cos(ra(0 - ip E )) 




where (• • •) denotes averaging over events. We define the second moment as 




and the fluctuations are 




There are several sources of trivial Vi fluctuations which are not related to the system's 
dynamics of interest. I start with the fluctuations caused by a varying number of particles used 
to determine tpE an d 1)%. I assume here that the Fourier amplitudes v n do not change from 
event to event and that the only correlations in the system are those due to the flow. Then, the 
azimuthal distribution of N particles is a product of iV single particle distributions 

P£((t>l,(j>2,---,<t>N)=V N P ev (<t>l)P ev ((f>2) ■■■ PeAM, (11) 

where Vn is the multiplicity distribution while all distributions P ev (4>i) are given by Eq. pi)|). 
The single particle distributions P ev (4>i) are correlated to each other because of the common 
angle ip R . 

Using the distribution (|11|) . one finds (neglecting ^4) the variance of V2 as 

'2\ _ m \2 



1 , , v2 (R 2 2 ) - (R 



Var ^ = ww +( ^' 1 w^ (12) 

which holds for (N) 3> 1 and small multiplicity fluctuations. The second term in r.h.s of Eq. (|12j) 
appears to be much smaller than the first one as {R%) — (R2} 2 ~ {M)~ 2 , where M is the number 
of particles used to determine the reaction plane and M is assumed to be of the same order as 
N. Thus, the statistical noise contribution to 5v 2 = ^/Var(t>2) is finally estimated as 

* 2 - jmm ■ (13) 

As well known, (i^) strongly depends on the collision impact parameter b. Using the 
parameterisation of this dependence given in 0, one computes 5v2 as (d(v 2 ) /db)5b. The impact 
parameter is measurable through the multiplicity N p of participating nucleons. N p is directly 
related to b. For b ~ 10 fm, when the flow in Au-Au collisions is maximal, the V2 fluctuations 
due to the impact parameter variation vanish because d(v2)/db = 0. For 6 = 5 fm, where 
(v 2 ) ~ 0.03, one finds 5v 2 ~ 8 • 10~ 4 5N p . When SN P = 30 and (N) = 500, the magnitude of the 
V2 fluctuations caused by the impact parameter variation is approximately equal to that of the 
statistical noise. 

The next source of trivial V2 fluctuations is a variation of thermodynamic parameters. The 
contribution caused by the particle number fluctuations can be estimated as 

where the effective power P is 

_ _ dln(v 2 ) (N) d(v 2 ) 



dln(N) (v 2 ) d{N) ' 
Assuming the poissonian character of multiplicity fluctuations, one obtains 

The value of the index P can be estimated within the hydrodynamic model j^Hj as P ~ 0.4. 
Comparing Eqs. ()13() and (|14[). one finds that the ratio of the thermodynamic fluctuations to 
the statistical noise, is 0.04 for (i^) = 0.07 and P = 0.4. Thus, the thermodynamic fluctuations 
are much smaller than the statistical noise. 

Concluding this section, I propose to perform a systematic measurement of event-by-event V2 
fluctuations for the centrality corresponding to the maximal elliptic flow. Then, the fluctuations 
caused by the impact parameter variation vanish or at least they are very small. If the flow 
is built up in the course of hydrodynamic evolution of the equilibrium system, 1)2 should be 
dominated by the statistical noise related to the finite particle number. The noise can be 
identified due to the characteristic 1/{N) dependence. 



3. 2. $ measure of azimuthal fluctuations 

Since a measurement of v-i fluctuations discussed in the previous section is rather difficult, I 
suggest to consider a much simpler integral measurement of azimuthal fluctuations [H2], using 
the so-called $ measure introduced in |35| . 

The correlation (or fluctuation) measure $ is defined as follows. One defines the variable 

z = f x — x, where x is a single particle characteristics such as the particle transverse momentum 
or the azimuthal angle. In this section the overline does not denote averaging over particles 
from a single event but averaging over a single particle inclusive distribution, x is identified here 
with the particle azimuthal angle. The event variable Z, which is a multiparticle analog of z, is 
defined as Z = f J2iLi( x i ~ x), where the summation runs over particles from a given event. By 
construction, (Z) = 0. The measure $ is finally defined as 



^ def 



(N) 



$ obviously vanishes in the absence of any inter-particle correlations. Other properties of are 
discussed in [Hfij . 

The $ measure is sensitive to the azimthal fluctuations caused by the transverse collective 
flow. Let me compute it, assuming that the only correlations present in the system are due to 
the collective flow. The inclusive 4> distrubtion, which is flat in the range [0, 2tt], provides 4> = tt 
and 4> 2 = Tin 2 , and consequently, z 2 = i-7r 2 . Since Z = X)i=i(^i ~~ 0)> one computes (Z 2 ) 



the event distribution (|lUj) . as 

r2n i-2-k r 2ir 



using 



(Z 2 ) = -^J2V N #!.../ d^ N P^fa) . . . P ev (07v) (01 + • • • + <PN ~ Ncf)) 2 . (15) 

JO 27T ^ Jo Jo 



The averaging over the amplitudes v n , which is not shown here, is implied. At first glance, 
the multi-particle distribution from Eq. (|15j) might look as a simple product of one-particle 
distributions. One should note however that every P e v(^) depends on the reaction plane angle 
ipfi, and the integration over ipR leads to the correlated multi-particle distribution. 

After elementary calculation, one finds for small v n , poissonian multiplicity distribution and 
(N) ^> 1, an approximate expression of interest 
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If all amplitudes v n except V2 vanish, as it approximately happens in the central rapidity domain, 
and V2 equals a unique value 0.07, Eq. ((ToT) for (N) = 500 gives $ = 0.2. 

As already mentioned, the transverse flow is far not the only source of the azimuthal 
fluctuations. In particular, the effect of quantum statistics contribute here. To estimate the 
effect one computes $ in the equilibrium ideal quantum gas. The result reads |32) 
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A denotes the fugacity, the upper sign is for fermions and the lower one for bosons. As seen, <I> is 
independent of the system's volume and of the number of particle's internal degrees of freedom. 



For massless bosons with vanishing chemical potential, Eq. (|17j) gives ~ 0.3 for any T. More 
realistic calculations |H2j provide ~ 0.06 for chemically equilibrated pions at T = 150 MeV. 

I conclude this section by saying that a measured value of <J>, which would significantly exceed 
predictions of Eq. (|16|) with non- fluctuating amplitudes v n , would be an obvious signal of sizable 
dynamical fluctuations. 



4. Final remarks 

The magnetic instabilities provide a plausible mechanism responsible for a surprisingly short 
equilibration time observed in relativistic heavy- ion collisions. Fast isotropization is a distinctive 
feature of the mechanism. It is certainly desirable to look for experimentally detectable signals 
of the instabilities driven thermalization. In my talk I have proposed to study azimuthal 
fluctuations, in particular the event-by-event fluctuations of the elliptic flow which is generated 
at the collision early stage. I have not been able to present a quantitative prediction but 
observation of sizeable dynamical fluctuations would be a strong argument that behind a smooth 
hydro dynamic evolution there is a violent phenomenon of plasma instabilities. 
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